Disorder and relaxation mode in the lattice dynamics of 
PbMgi/3Nb2/303 relaxor ferroelectric 
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The low-energy part of vibration spectrum in PbMgi/3Nb2/303 relaxor ferroelectric was studied 
by inelastic neutron scattering. We observed the coexistence of a resolution-limited central peak 
with strong quasielastic scattering. The line- width of the quasielastic component follows a Fq -I- Dq^ 
dependence. We find that Fq is temperature-dependent. The relaxation time follows the Arrhenius 
law well. The presence of a relaxation mode associated with quasi-elastic scattering in PMN indicates 
that order-disorder behaviour plays an important role in the dynamics of difi^use phase transitions. 

PACS numbers: 77.80.-e, 61.12.-q, 63.50.4-x, 64.60.-i 
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Enormous attention has been paid to relaxor ferro- 
electrics for more than 40 years due to their intriguing 
physical properties. A marking feature observed in re- 
laxors is the frequency-dependent peak in the dielectric 
permeability e which typically extends over a few 100 
degrees and does not link directly to the macroscopic 
changes of symmetry Q . Since many physical properties 
of relaxors exhibit anomalies in this temperature range, it 
was called " diffuse phase transition" . Despite intense 
investigations, however, the nature of the diffuse phase 
transition in these compounds remains unclear ^|^. 

The well-known crystal PbMgi/3Nb2/3 03 (PMN) is a 
model relaxor. In PMN, the peak in the real part of the 
dielectric permeability e' appears around the mean Curie 
temperature Tcm ^ 270 K |3] . By application of an elec- 
tric field it is possible to induce a structural phase tran- 
sition at T~ 210 K 0. Furthermore, at Td ~ 620 K the 
temperature dependence of the optical refraction index 
deviates from the expected linear dependence 0. This 
result was explained by the appearance, far above T,„ , of 
small polarized regions within the crystal . Thispolar- 
ization would exist on the scale of a few unit cells |3l only, 
and hence called 'polar nanoregions' (PNR). The concept 
of PNR became very popular to account for properties 
like the diffuse scattering observed by neutron scattering 
in relaxors of PbB'^/gB^'/gOa type |ili,ESIi3- 

To obtain a theoretical description of a phase transi- 
tion, an important point is to know whether this tran- 
sition is primarily of displacement, or of order-disorder 
type Since PMN crystallizes in the perovskite 

structure where the B-site is disordered, it was expected 
that the diffuse phase transition is of displacement type. 
Within such a scenario a soft-phonon mode should be 
observed. Extensive search for the presence of a soft 
mode in PMN was performed both by light (for a review 



see s-nd neutron scattering. From light scattering no 
soft mode has ever been found. On the other hand, neu- 
tron scattering experiments revealed a very complicated 
behaviour of the phonons in PMN [ij, \\A I1(t| . Firstly, 
a soft excitation was identified in PMN above T^ [lj|, 
which is coupled with transverse acoustic phonon. Later, 
in Ref. |l5l | a strong renormalization of the lowest trans- 
verse optic (TO) phonon branch was observed start- 
ing from T~1100 K. Recently, a soft mode coupled to 
the lowest TO phonon was reported Obviously, 
the identification of a soft mode in PMN has not been 
definitely settled. In addition to the phonon modes a 
resolution-limited central peak (CP) is observed by neu- 
tron scattering in PMN without and with an applied 
electric field [l3|. The timescale of CP is also at the 
center of a controversy because it contradicts light scat- 
tering data [l^ . where a broad quasielastic component 
was observed with a typical width of 2 meV. We note, 
that the search for a broad quasielastic component in re- 
lated relaxor PbZni/3Nb2/303+8%PbTi03 was not suc- 
cessful jl9| which again is in contradiction with results 
of light scattering {e.g. 113). Hence, neutron scattering 
data speaks in favor of displacement behaviour in relaxor 
ferroelectrics coexisting with a narrow resolution limited 
CP. On the contrary, light scattering data suggest order- 
disorder behaviour in these compounds. 

To solve these discrepancies between light and neu- 
tron scattering results and to show that order-disorder 
behaviour plays an important role in the lattice dynam- 
ics of PMN we have studied the low-energy range of 
the vibration spectrum in this crystal by inelastic cold- 
neutron scattering under improved resolution conditions 
with respect to the previous studies of refs. 0,0,0]. 
The measurements were carried out with the three-axis 
spectrometer TASP, located at the neutron spallation 
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source SINQ at the Paul-Scherrer-Institutc, Switzerland. 
A high-quality single crystal of PMN 8 cm'^, mo- 
saic ^ 20') was mounted inside a closed-cycle refriger- 
ator equipped with a small furnace. The crystal was 
aligned in the (h h l)-scattering plane. The measure- 
ments were performed in the temperature range 100 K 
- 450 K. The (002) reflection of pyrolytic graphite (PG) 
was used to monochromate and analyze the incident and 
scattered neutron beams, respectively. The spectrome- 
ter was operated in the constant final-energy mode with 
either kf = 1.64 A"! or kf = 1.97 A'^ using a PG fil- 
ter to remove higher-order wavelengths. The horizontal 
coUimation was guide— 80' — 80' — 80'. With that config- 
uration the energy resolution as measured with a stan- 
dard Vanadium sample is 0.4 mcV for kf ~ 1.97 
and 0.2 meV for kf = 1.64 A~^. Figurenshows {q, huj) 
contour-maps of neutron scattering spectra in the vicin- 
ity of Bragg positions (0 1), (1 1 0), and (0 2) at 
T = 300 K, respectively. It is clear that all three fig- 
ures differ from each other. For example, the intensity 
distribution in the (0 1) Brillouin zone (BZ) consists 
of a Bragg peak, a rather intense signal which is narrow 
in energy and broad in g, and a diffuse signal which is 
broad both in q and in hu. In the (0,0,1) BZ no linear 
dispersion which could be associated with TA phonons 
is observed. On the other hand, one can easily sec the 
linear dispersion of TA phonons in the (0 2) BZ, in 
addition to the Bragg peak and to the intense signal nar- 
row in energy. This signal corresponds to the central 
peak reported before [l4 . . Qualitatively, the distri- 
bution of the intensity in the (1 1 0) BZ lies between 
what is observed in the (0 1) and (0 2) zones. The 
typical peak intensity of a TA phonon measured close the 
(0 2) Brillouin zone center is ~ 5000 counts at room 
temperature. It is well-known that the neutron intensity 
of an acoustic phonon scales with the structure factor of 
the corresponding Bragg peak and with the square of the 
modulus of the scattering vector . The structure 

factor of the (0 2) Bragg peak is 100 times larger then 
that of (0 1). Also, increases by a factor of ~ 4 while 
measuring in the (0 2) BZ. Thus, the total decrease of 
intensity of the TA phonons in (0 1) BZ with respect to 
(0 2) amounts to 400. Hence it is impossible to detect 
the intensity of such a phonon with respect to the diffuse 
intensity which amounts to ^^200 counts in the (0 1) 
zone. 

In this letter we will concentrate on the properties of 
the quasielastic intensity found in PMN and not on the 
behavior of the acoustic phonons as these have been the 
subject of numerous studies [1, 0, 0, 113 • In order 
to interpret the data quantitatively we decomposed the 
neutron scattering spectra /(Q,cj) in the following way: 

/(Q, u) - [Sep + Sq^ei. + Sdho] ® R{Q: u;) + B (1) 

where Sep refers to the resolution- limited central peak; 
Sq-ei. is the quasi-elastic scattering shown in Fig. 1; 



0.15 
Did 
OlOJ 
-OJW 
■OiUS 
-O.Jfl 









B.10 
















-O.Ofl 


1 — 1 










-ans 








0.\5 




Mfl 










£1. 






■OJM 






a; 










-O.lfl 




-IM -l.OO tm IJM 2,00 3.00 

Energy (meV) 



FIG. 1: Color maps of neutron scattering intensity around 
the (001) (top), (110) {middle), and (002) {bottom) Bragg po- 
sitions at T = 300 K. Asymetry in the intensity contours of 
neutron energy gain (negative energies) and neutron energy 
loss (positive energies) parts of spectra is due to the resolu- 
tion effects. Note that the intensity is given in a logarithmic 
scale. 1 rlu=1.555 . 



Sdho describes the phonon scattering. The symbol ® 
stands for the convolution with the spectrometer resolu- 
tion function i?(Q,a;) B denotes the background 
level. 

We found that the data could be modeled using 



ScP^A{q)5{u), 



Sq-el = I{q) 



1 r„ 



l-exp(-tj/r)7rw2 + r2' 



(2) 
(3) 
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and 

SDHO = W{q)- JTf^-r^ f^]s2 2~2- 

Here 5{uj) is the Dirac delta function; Vq and 7^ are the 
half-width at half maximum (HWHM) of the Lorentzian 
and damped-harmonic oscillator function, respectively. 
A{q), I{q) are the q-dependent susceptibilities of the cor- 
responding excitations and W{q) is the TA dynamical 
structure factor. The renormalized dispersion of the TA 

phonons is given by Vlq = yG^^Tl^ where ujq is the fre- 
quency of a damped oscillation. In our description, f2g 
is the physically relevant quantity |24| . Since the mea- 
surements are restricted to small values of q, we used a 
linear dispersion for the TA phonon. The validity of such 
an approximation is clear from the intensity distribution 
shown in Fig. 1. Fig. [21 shows typical constant-Q scans 
in the (0 1) and (1 1 0) BZ, respectively. Comparison 
between Figs.Et audits shows that phonons are visible 
at 7ia> = ± 1.6 meV for Q=(l 1 0.075). On the contrary 
there is no phonon scattering at Q=(0.05 0.05 1). Hence, 
to analyze the data in the (0 1) BZ, we fixed W{q) to 
zero. The solid lines in Fig.|21show that the data are well 
represented by equation ^ Figs.|21i,b show the q- depcn- 




^ 2 4 6 8 10 12 




2 4 6 8 10 12 
Energy (meV) 

FIG. 2: Typical examples of constant-Q scans in different 
BZ. The bold line shows the results of data analysis as de- 
scribed in the text. The solid line shows the contribution of 
the quasielastic component, the dash-dotted line the narrow 
central peak and the dotted line the TA phonon. 

dence of the integrated intensity and of the damping of 
the quasielastic scattering observed in the (0 1) BZ at 
T = 450 K and 225 K. The quasielastic scattering can 
be approximated by a Lorentzian line-shape in q-space. 
Then the correlation length is ^ = 6 A at T = 450 K, and 
^ = 12 A at T = 225 K (^=1/HWHM of the Lorentzian). 
The damping of the QE component is found to behave 



like V{q) = Fq + Dq^ . At T = 450 K the parameters 
are Fq = 0.28 ± 0.01 meV and D = 24.5 ± 7.3 meV-A^ 
and at T = 225 K we obtain Tq = 0.02 ± 0.01 meV and 
D = 21.3 ± 8.9 meV-A^. The absence at all q of an in- 
trinsic energy width of CP implies that the associated 
fluctuations are very slow. Further, the intensity of the 
QE scattering does not scale with the narrow CP. In the 
(0 2) BZ (FiglU bottom) the CP is very intense, but 
we were unable to detect QE component in this BZ. On 
the contrary, in the (0 1) BZ the narrow CP is signifi- 
cantly weaker and the QE component is clearly seen. The 
different q-dependences of the energy widths and behav- 
ior of the relative intensities of the CP and quasi-elastic 
scattering indicate that these processes are of different 
origins. 

We now turn to the temperature behaviour of the re- 
laxation time r = 1/F of the QE-mode measured in the 
(0,0,1) and (1,1,0) BZ, respectively, and shown in Fig.|3t. 
It is obvious, that the temperature dependence of r is 
very similar in both zones. We found that r is well de- 
scribed by the Arrhenius law t{T) = tq • exp [Ea/ks ■ T]. 
The solid line in Fig. |3Jd corresponds to the parameters 
To = 8.6 ± 0.1 • 10-12 and Ea = 0.023 ± 0.0002 (eV). 
While these parameters are reasonable, they have to be 
taken with care as quantities like the refraction index Ml, 
the specific heat [23 , and the diffuse scattering |^ |2g 
exhibit anomalous behaviours in this temperature range. 
This implies that the structure of PMN undergoes sig- 
nificant changes on a local scale, which in turn should 
affect the value of Ea ■ Despite this uncertainty, the fact 
remains that t{T) is connected to a thermally activated 
behaviour, which allows us to speculate about the origin 
of the relaxation mode. It is well-known, that in relax- 
ors, Pb ions are shifted from their ideal positions in the 
perovskite structure [2^. Recently it was shown that 
the probability density function p{R) for Pb in PMN has 
more than a single maximum below T^ |27j . Therefore, 
it is reasonable to assume that at high temperature the 
Pb ions move between maxima of p{R). In that case 
when the temperature is low enough with respect with 
the potential barrier, the hopping frequency l/r becomes 
small. This gives rise to a narrowing in the energy-width 
of the QE component, which, within the limitations of the 
present experiment, is resolution limited at T = 100 K 
in PMN. 

There is a number of models which describe the pres- 
ence of a central peak in the vicinity of a structural phase 
transition (for a review of earlier papers see [l^). How- 
ever, PMN does not undergo any structural phase tran- 
sition without application of an electric field. To link the 
coexistence of a central peak with quasielastic scattering 
in a broad temperature range might be of importance to 
understand the diffuse phase transition in relaxor fcrro- 
electrics. 

To summarize, we have carried out high-resolution 
neutron scattering measurements in the relaxor fcrroelec- 
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FIG. 3: a and b q- dependence of integrated intensity 
(squares) and of damping (triangles) of the QE component at 
T = 450 K and 225 K, respectively . c Temperature depen- 
dence of the relaxation time r associated with the quasielastic 
component. The solid line is a fit to the data with the Arrhe- 
nius law. 



trie PMN. The low-energy spectra exhibit two compo- 
nents associated with two different timescales. One is the 
resolution-limited central peak. The other component is 
quasi-elastic and appears to be related to ionic motion 
of thermally activated character. The damping of this 
relaxation mode follows a dependence. The presence 
of a relaxation mode associated with quasi-elastic scat- 
tering in PMN indicates that order-disorder behaviour 
plays an important role in the dynamics of diffuse phase 



transitions. It will be interesting to see if this relaxation 
mode, which was also observed in PMT , is a common 
feature of all relaxor ferroelectrics. 
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